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We present a new method of bosonization of fermion systems applicable when the partition func- 
tion is dominated by composite bosons. Restricting the partition function to such states we get an 
euclidean bosonic action from which we derive the Hamiltonian. Such a procedure respects all the 
fermion symmetries, in particular fermion number conservation, and provides a boson mapping of 
all fermion operators. 
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I. INTRODUCTION 

The importance of effective bosons in nuclear physics 
became clear after the observation that heavy deformed 
nuclei display some of the features of superconducting 
systems jjj. In these nuclei there must then be Cooper 
pairs of nucleons. In this spirit Arima and Iachello in- 
troduced two different composite bosons, the s- and 
d-bosons. Their model, the Interacting Boson Model, 
proved extremely successful in reproducing low energy 
nuclear properties, but has not yet been derived in a fully 
satisfactorily way from a nuclear Hamiltonian. 

Many attempts to reformulate the nuclear Hamiltonian 
in terms of effective bosons, however, were done before 
the Interacting Boson Model was invented. Of special 
importance are the works of Beliaev and Zelevinsky Q , 
who constructed a composite boson operator requiring 
its commutation with the nuclear Hamiltonian, and of 
Marumori, Yamamura and Togunaga Q, who developed 
a method based on a map of fermion into boson matrix 
elements. 

The first important step in the derivation of the IBM 
respecting nucleon number conservation is due to Otsuka, 
Arima and Iachello Q • Their work is based on a map of 
a single j-shell nucleon space into a boson space. The 
boson Hamiltonian so obtained reproduces exactly the 
spectrum of the pairing model. Their procedure has been 
somewhat extended |6j but not sufficiently generalized. 

There are several recipes for bosonization based on 
a mapping of the nucleon model space into a boson space. 
Such methods do not violate nucleon conservation and in 
principle yield an exact solution to the problem, but in 
practice one has to perform a truncation in the nucleon 
space related to a selection of degrees of freedom guided 
by physical insight and calculational convenience. One 
shortcoming of this procedure is the appearance of "in- 
truders" , namely states which in spite of their low energy 
do not appear in the boson space generated by the map- 
ping yj. 

A different approach to bosonization which does not 
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require a preliminary truncation of the nucleon space 
and does not violate nucleon conservation j^, ^| is based 
on the Hubbard-Stratonovich transformation. The latter 
renders quadratic the fermion interaction by introduc- 
ing bosonic auxiliary fields which in the end become the 
physical fields. The typical resulting structure is that of 
chiral theories In such an approach an energy scale 
emerges naturally, and only excitations of lower energy 
can be described by the auxiliary fields. In our opinion 
this approach has not received enough attention and its 
potentiality has not been fully explored. 

The physical idea behind bosonization is that certain 
composite bosons dominate the partition function at low 
energy (Boson Dominance), an assumption certainly jus- 
tified for Goldstone bosons like Cooper pairs. We present 
a new way to implement Boson Dominance. We intro- 
duce generic nucleon composites whose structure will be 
determined at the end by a variational procedure, and 
evaluate the partition function restricted to such com- 
posites. In this way we get an euclidean bosonic action 
in closed form. 

In the derivation of the effective action we need only 
one approximation, concerning the identity in the space 
of the composites, but we respect all the nucleon sym- 
metries, in particular nucleon number conservation. We 
emphasize that the closed form of the action opens the 
way to numerical simulations of fermion systems in terms 
of bosonic variables, avoiding the "sign problem" [l2| . 

Bosonization is achieved within the path integral for- 
malism. In this framework the standard procedure to 
evaluate physical quantities is to first find the minimum 
of the action at constant fields. Depending on the solu- 
tion, one has spherical or deformed nuclei. In the latter 
case rotational excitations appear as Goldstone modes as- 
sociated to the spontaneous breaking of rotational sym- 
metry. The notion of spontaneous symmetry breaking 
survives in fact with a precise definition also in finite 
systems |lCj. Next the quantum fluctuations must be 
taken into account. But we anticipate a subtlety in this 
program which in the present case is not completely stan- 
dard due to the composite nature of the bosons. 

In nuclear physics the Hamiltonian formalism is of 
much wider application. Since the effective bosonic ac- 
tion, due to compositness, is not in canonical form, it 
has been necessary to devise an appropriate procedure 
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to derive the Hamiltonian. In this context the mentioned 
subtlety finds a natural solution. 

Only composites which have components on many nu- 
cleon states can be approximated by bosons. But the bo- 
son space cannot be arbitrarily truncated. For instance, 
even if the nuclear potential contains only monopole and 
quadrupole pairing interactions, the s- and d-bosons will 
be coupled to all the other bosons permitted by angu- 
lar momentum conservation. An approximate decoupling 
should arise dynamically, but at the moment a clearcut 
mechanism is not known. We will say more about this 
point in our Conclusions. 

Bosonization appears in several many-fermion systems 
and relativistic field theories. The effective bosons fall 
into two categories, depending on their fermion number. 
The Cooper pairs of the BCS model of superconductiv- 
ity, of the Interacting Boson Model of Nuclear Physics, of 
the Hubbard model of high T c superconductivity and 
of color superconductivity in QCD have fermion num- 
ber 2. Similar composite bosons with fermion number 
zero appear as phonons, spin waves and chiral mesons in 
QCD. The latter bosons can be included in the present 
formalism by replacing in the composites one fermion 
operator by an antifermion (hole) one. This becomes 
necessary when the interaction contains, as in nuclear 
physics, important particle-hole terms. This extension of 
the method will be presented in a future work. 



II. OUTLINE OF THE APPROACH 



Consider a nuclear partition function 



Z = tr 



> x|> j --^(H - fi N n N ) 



(1) 



where T is the temperature, [An the nucleon chemical 
potential and the nucleon number operator. A sector 
of tin nucleons can be selected by the constraint 



T 



d 



■\uZ = TIN- 



N 



(2) 



Under the assumption of Boson Dominance we can re- 
strict the trace to nucleon bosonic composites. The re- 
stricted partition function can be written 



tr 



Pcxp ( -7j;{H - HNn N ) 



(3) 



where V is a projection operator in the subspace of the 
composites. We will only be able to construct an ap- 
proximation to such operator. This is the only approx- 
imation we will do (beyond the physical assumption of 
Boson Dominance). In analogy to the case of elementary 
bosons we will assume 



V= / dn(f3*,(3)\(3)(/3\ 



(4) 



A preliminary presentation of our results was given in 
0. 

The paper is organized as follows. In Section II we out- 
line our approach. In Section III we report the bosonic 
effective action. In Section IV we derive the bosonic 
Hamiltonian as the normal ordered form of a nonpolyno- 
mial function of creation-annihilation boson operators. 
This result includes the boson mapping of all fermion 
operators. For a practical use this Hamiltonian must be 
expanded in the inverse of the dimension of the nucleon 
space. In Section V we report an independent deriva- 
tion of the boson Hamiltonian valid for a small number 
of nucleons and we specify this Hamiltonian to the case 
a single j-shell. In Section VI come back to the path 
integral formalism, introducing the Goldstone and Higgs 
fields, and in Section VII we end with our conclusions. 

In the presentation of our results, to facilitate the un- 
derstanding of the logical development, we choose to 
relegate many technical details in a number of Appen- 
dices. In Appendix lAl we reported some results concern- 
ing Berezin integrals and in Appendix iBl their use in cal- 
culations with coherent states of composites. The basics 
of this formalism can be found in a condensed form in 
|15| , while an exhaustive presentation is given in |16| . In 
Appendix [C] we discussed the properties of an operator 
which approximates the identity in the space of the com- 
posites and in Appendix [D] some intermediate steps in 
the derivation of the effective bosonic action. 



where d/i(/3*, (3) is an integration measure (to be specified 
later) over the holomorfic variables /?*,/? and the \f3) are 
coherent states of composites 



|exp 



(5) 



Definition and properties of the operator V are discussed 
in Appendixlcl The \(3) are defined in terms of composite 
creation operators 



i 



2VTl, 



-J By = 



i 
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(6) 

The c''s are nucleon creation operators, the m's represent 
all the nucleon quantum numbers, the matrices Bj are 
the form factors of the composites with quantum num- 
bers J. O j is the index of nilpotency of the J-compositc, 
which is defined as the largest integer such that 



bj) ¥=0. 



(7) 



In the present paper we will assume for simplicity the 
index of nilpotency independent of the quantum num- 
bers of the composites, and equal to half the dimension 
of the nucleon space, but we will mention a possible con- 
sequence of this simplification. 

It is obvious that a necessary condition for a compos- 
ite to resemble an elementary boson, is that its index of 
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nilpotency be large. But such condition is in general not 
sufficient. Consider for instance the case 



= 6. 



m ^mi i 



(8) 



where d\ = 1, d m l,m ^ 1. Such a composite, irre- 
spective of its index of nilpotency, consists essentially of 
a unique state of a nucleon pair. We must instead re- 
quire that the composites actually live in a large part of 
the nucleon space. This can be ensured by the further 
requirement 



det 



1. 



(9) 



Solutions to the equations for the £?-matrices which do 
not satisfy the above condition must be discarded. 
Evaluation of the trace (reported in Appendix[D| gives 



d(3*d(3 

27TJ 



exp i 



-S eS (p*,(3)). (10) 



Bosonization is thus achieved, and the nuclear dynam- 
ics can be studied by functional or numerical methods. 
The last possibility appears interesting because it avoids 
the "sign" problem [l^] which affects the Monte Carlo 
approach to the study of many-fermion systems.. 

In nuclear physics it is generally used a Hamiltonian 
formalism. The Hamiltonian of the effective bosons, Hb, 
cannot be read directly from the effective action, because 
S e Q , (3) does not have the form of an action of elemen- 
tary bosons. Indeed it contains anomalous time deriva- 
tive terms, anomalous couplings of the chemical potential 
and nonpolynomial interactions, which are all features of 
compositeness. Therefore it has been necessary to de- 
vise an appropriate procedure to derive Hb 1 which also 
has been given in closed form, in terms of boson opera- 
tors b\b, satisfying canonical commutation relations. In 
conclusion 



(11) 



where \ib is the boson chemical potential and h B the 
boson number operator. 

For a practical use, however, it is necessary to perform 
an expansion of Hb- The expansion parameter is the 
inverse of fi. 



III. THE EFFECTIVE BOSONIC ACTION 

We write the nucleon-nucleon potential as a sum of 
multipole pairing terms, so that the Hamiltonian has the 
form 



H = ^h c-Y^9 K \^F^\cF KI 



(12) 



K 



The one body term includes the single-particle energy 
with matrix e, the nucleon chemical potential /ijv and any 



single particle interaction with external fields included in 
the matrix M. 



fi N + M. 



(13) 



The matrices Fk are the form factors of the potential, 
normalized according to 



ix{Fj Ci F Ka ) = 2SlS Kl K a 



(14) 



Any potential can be written in the above form [13 . But 
this form is not convenient, as it is well known, when 
particle- hole terms are important. To properly account 
for such terms in the present scheme it is necessary to 
introduce phonons, which will be done in a separate work. 

In order to evaluate Zq we divide the inverse temper- 
ature in Ao intervals of size r 



T 



Not 



(15) 



Then as shown in Appendix D the cuclidean effective 
action has the form 



S eS (f3* , 0) = N In J + t J2 tr j ^ In [1 + r r V, 
~\ X>* W^Fk) tr(ri^$) + 2 ((T - 1) F^Fk 



4>] 



A" 



-\r$*Fk,TF K $\. 



[r$* (<z>h T + h<s>) 



(16) 



where J is a function appearing in the measure defining 
the operator V, 



K 



V t / = - (ft+i - ft) , 
r 



(17) 



(18) 



(19) 



(20) 



and [..,..]+ is an anticommutator. Notice in the second 
line a trace inside the trace. The variables /?*,/? are al- 
ways understood at times t, t — 1 respectively. S e g has a 
global U(l) symmetry which implies boson conservation. 

The fermion interactions with external fields are ex- 
pressed in terms of the bosonic terms which involve the 
matrix M. (appearing in h) . 

The dynamical problem of the interacting (composite) 
bosons can be solved within the path integral formalism. 

Part of the dynamical problem is the determination 
of the structure matrices Bj. This can be done by ex- 
pressing the energies in terms of the Bj and applying a 
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variational procedure which gives rise to an eigenvalue 
equation. 

S e Q must be compared to the action of elementary 
bosons. If Hsib^b) is the Hamiltonian of these bosons 
in normal form, the corresponding action is |l5j 

Sb =TY,{(3*^t(3-H(P*,P)+n B (3*P} (21) 



where again the variables (3* , (3 are understood at times 
t, t — 1 respectively. 

We notice that S e g differs from Sb in many respects 

i) there is no canonical time derivative term 

ii) the coupling of the chemical potential ( appearing 
in h is also noncanonical ) 

hi) there are non polynomial interactions because of 
the T-function. This function becomes singular, as it will 
become clear in the sequel, when the number of bosons 
is of order CI, reflecting the Pauli principle. 



only after an appropriate subtraction, which can be per- 
formed by the change of variables 



K • 



(25) 



where r is a parameter which will be fixed later. This 
subtraction corresponds to the Bogoliubov transformation 
in other approaches, but does not violate nucleon number 
conservation. We can then rewrite the T-function in the 
form 



r 2 r', 



where 



r' = 



i 



1 + - (/?'* B [3' -B ] - 



(26) 



(27) 



IV. THE BOSONIC HAMILTONIAN 

Let us start by examining the features of composit- 
ness when the number of bosons is much smaller than 
fl. Since the expectation value of (3* ■ (3 is of the order 
of the number of bosons, in this case we can perform an 
expansion of logarithm and T-function in inverse powers 
of Q. Expanding the logarithm we have 

Y tr In [1 + rr$* V$] = -tr ($* V t $) 



-tr [$*$$*V t $l 
4 



(22) 



The first term can be made canonical by normalizing the 
boson form factors as the potential form factors 



tY{B\B K ) = 2flS, 



J,K- 



(23) 



The other terms are then of order fi -1 . Notice that the 
diagonal condition is only a matter of normalization, but 
the off diagonal one must be compatible with the dy- 
namics. If this is not the case, a redefinition of the (3's is 
necessary. 

Expanding the T-function we get the following cou- 
plings of the nucleon chemical potential 



UN tr 



(24) 



Only the first term is canonical, and /ijv is not half the 
boson chemical potential as one might expect. As we 
will see in Sec. V, for n <C Q these anomalous couplings 
can be eliminated by a redefinition of the chemical po- 
tential, so that in the case of a small number of bosons 
the Hamiltonian can be derived without difficulty. 

But when the number of bosons is of order f2, the ex- 
pansion of logarithm and T-function can be performed 



For a suitable choice of r, T' admits an expansion in . 
Now we take the function J appearing in S b q equal to 
the jacobian of the transformation (|25|l 



-N B 



(28) 



where Nb is the number of bosonic degrees of freedom. 
Therefore the partition function becomes 



Zc = 



d(3'*d(3' 



2m 



exp(~S") 



(29) 



where 



S\(3'* , 13') = t tr | ^ ln C 1 + r V V* 



2 (T'F^Fk) - [T'&*F^ K ,T'F K &] 



\ [r'$'* ($' h T + h $')] +\Y,9k (fIf k ) \ 

K ) 



with 



(30) 



(31) 



We assume, and we will verify a posteriori, that the pa- 
rameter r can be chosen in such a way that the anomalous 
time derivative terms be of order f2 -1 . Then to this order 
Zq can be written as a trace in a boson space 



(32) 
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The Hamiltonian H' is obtained [l^ by omitting time 
derivative and chemical potential terms, and replac- 
ing the variables /?'* , (3' by corresponding creation- 
annihilation operators , b. These satisfy canonical com- 
mutation relations and should not be confused with the 
corresponding operators for the composites which are dis- 
tinguished by a hat 



H'(r, m ) = : tr J -j£,9k (l ~ 77 



K 



x tv{T b F K $ b ) + 2 T b F^Fk - [T b ^ K> T b F K ^ b l 

+~ [r b ^ ($ 6 ^ T + h $„)] + ~ J] 0x4^* I : ■ 



K 



The colons denote normal ordering and 



(33) 



1 + - (V ■ B 6 • 5+ 



(34) 



Here we meet with a subtlety. H' commutes with the 
boson number operator, so we can select sectors with 
a given number of bosons. But we are not guaranteed 
that these bosons carry nucleon number 2, because of 
the noncanonical coupling of the chemical potential. We 
can enforce this fundamental property in the following 
way. Let us denote by E' (n) the lowest eigenvalue of H' 
in the sector of n bosons. We require that E' (n) be the 
lowest eigenvalue for n = 5 njv 



d 1 
E' \ = 0, n = -n N . 



dri 



(35) 



This determines r as a function of the number of bosons 
and the nucleon chemical potential: r = r(n, fAjsr), ensur- 
ing that the bosons carry nucleon number 2. Condition 
P| then determines the nucleon chemical potential as a 
function of n: Ji N = Jl N (n). The boson Hamiltonian in 
the sector of n bosons is finally 



H B {n) 



H'(r,n N ) 



2/j, N n. 



(36) 



It depends on n explicitly and through the dependence 
on n of r, ~p N . Therefore also the matrices Bj will depend 
on n, namely the form factors of the bosons depend on 
the number of the nucleons. 

Notice that H 1 provides the mapping of the nucleon 
interactions with external fields 



1 



JMc^: -tr [T b ^* b (^ b M T +M * 6 )] 



(37) 



s-boson condensates 



This procedure becomes particularly simple if the 
ground state contains only one species of bosons, which, 



having in mind the Interacting Boson Model, we call s- 
bosons. We will refer to such a ground state as to an 
s-boson condensate". In such a case, assuming for sim- 
plicity 



F$F = BlBo = 1, 



(38) 



in the evaluation of the ground state energy we can set 
1 



$b = 

r b = 



boBl 



1 



-1 -1 



(39) 



It is then convenient to adopt the following definitions 
1 



9k 



4fi 2 



|tr(i?oi4)|V 



9k 



K 



G = Y.9K- 

K 



(40) 



Let us disregard for a moment normal ordering. This is 
equivalent to the semiclassical approximation in the path 
integral formalism, namely to neglect quantum fluctua- 
tions. We then get for the lowest eigenvalue of H' in the 
n boson sector 



l + I(n-r 2 ) 



QG 1 



r 

IT 



2 hn 



-QG- ( l + -(n-r 2 ) 



1 — — J n(QG-G 



(41) 



where 



h = — tr h. (42) 
Condition i|35[) determines r as a function of njv, 

_ 2 _ -2(fi + n)h + fi(0 - 2n)G - 20G 
Q(-2h- 2G + flG) 

Inserting this value in E' a we get 



K = - 



riG-G 

Condition J5J then determines /in 
1 



h- -(QG-2G) 



Mat 



where 



e- -CiG + nG- -G 



(44) 



(45) 



(46) 
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It is not surprising that with this value of /i N 



r 2 = n. 



(47) 



We then get also the lowest eigenvalue of Hb (neglecting 
normal ordering) 

E = E' Q + 2n N n = 2ne-nnG + n 2 (g - i G^j . (48) 

In the case of a monopole pairing interaction, G = G = 
go, comparing to the exact spectrum (|65|l . we see that 
the coefficients of the powers of n are affected by errors of 
order f2 _1 . The form factor of the s-boson is determined 
by the minimizing Eq 



Bn — Fn. 



(49) 



In order to get an expression of Hb of practical use, we 
must perform an expansion in fl^ 1 . Since the energy 
scale is set by the single-particle energies, we must make 
an assumption concerning the magnitude of the coupling 
constants gx with respect to e. For a system with in- 
finitely many degrees of freedom, — > oo, in order to get 
finite energies we must assume gx ~ _1 , in which case 
/ijv is of order DP. Such a behavior is also acceptable for 
many nuclei. 

To take into account quantum effects we must put H' 
in normal order. This corresponds to include quantum 
fluctuations in the path integral formalism, and requires 
an expansion with respect to f2 _1 . Remaining in the case 
in which the ground state is an s-boson condensate, we 
only need the following equation 

1 oo 

-<6 " ]>>:(6t6 ) S : (&$)") 

n 

= ^2c s n(n -l)...(n- s + 1). (50) 

s=0 

Notice that the expectation value of any normal ordered 
function in a state of n bosons is a polynomial of degree 
not greater than n. 

A further simplification occurs if the number of the 
other bosons, which we denote by the label K, is much 
smaller than Q. In such a case we can obviously classify 
the terms appearing in the function r^, Eq. (|34f> . accord- 
ing to 



btb -r 2 +J2 b^B^B^ ~ 1 

Kx,K 2 

blBoJ^hc^ ~ ( 51 ) 



K 



It is then easy to see that neglecting terms of order 
or smaller, Hb is at most quartic in the K boson opera- 
tors. 



V. AN ALTERNATIVE DERIVATION OF THE 
HAMILTONIAN FOR n < fi 

In this Section we restrict ourselves to the case of a 
small number of bosons. Then the subtraction is not 
necessary, we can set r = 0, and we can put the effec- 
tive action in canonical form by a shift of the chemical 
potential. It is to emphasized that no other quantum 
corrections are necessary after such a shift. 

For simplicity we assume the coupling go to be positive 
(attractive pairing force) and larger than the other ones, 
so that at the minimum only /3o is different from zero. 
Assuming Bo to satisfy Ea. (|38|l . S e g- at constant fields is 

TS = fi^ 5i ,-r!^ |/? | 2 + ^.g^ ( 1 + 7^°l 



+^tr(e-Miv)|/3o| 2 ( 1 + ^\Po\ 2 



(52) 



Its minimum with respect to |/?o| 2 can be determined ex- 
actly, but since we will perform the I/O expansion we put 
ourselves in this framework since the beginning. We will 
retain only the first order corrections, which are of or- 
der 17°, with the exception of the coupling with external 
fields which are of order O -1 . We remind that the diffi- 
culties in the derivation of the boson Hamiltonian are due 
to anomalous time derivative terms and couplings of the 
chemical potential. In this approximation the first diffi- 
culty is overcome because, as already noted, noncanonical 
time derivatives are of order I/O. It remains to get rid 
of the noncanonical couplings of the chemical potential. 
For this purpose we set 



fi N = f2/ii + -hb 



(53) 



and expand wr to 1 /Q 



TS = -0(2 A11 + 5o )|/3o| 2 + (2e- AlB )|/3 | 2 + 2( Ail + 5o )|/9o| 4 . 

(54) 

Since \(3o\ 2 ^ 0, at the minimum the first term must 
vanish separately from the others and we get 



1 



IA,| 2 = -i\»B-e) 

go 2 



(55) 



We select a sector with a given number n of bosons by 
imposing the condition 



d 
dfj, B 



(TS) = n 



which yields 



|/?o| 2 = n. 



(56) 



(57) 



We see from Eq. Q24H that to order 0° the only noncanon- 
ical term is proportional to /Ui, which does not depend 



7 



on n. We can then insert in the action the definition l|53|) 
and get a canonical bosonic action with canonical chem- 
ical potential [Ib- 

There remains a last point. The energies are given by 
(minus) the logarithm of the partition function plus the 
chemical potential times the number of bosons. But we 
can subtract from the action the term f2/ii/3j • 0t—i, and 
subtract in the end from the energy only [Ib times the 
number of bosons. In this way we get exactly the boson 
Hamiltonian Hb- 



A. Few bosons in a single j-shell 

If the nucleons live in a single j-shell the form factors 
of the composites either vanish, or are equal to the form 
factors of the potential. But unless I << j condition @ 
is not satisfied, the composites do not have a high index 
of nilpotency and must be excluded. 

We identify the quantum number K with the boson 
angular momentum, K = (Ik, Mk), so that the form fac- 
tors of the potential are proportional to Clebsh-Gordan 
coefficients 

{F m )m u ma = V2Q (j mi jm 2 \IM) , Q = j + -, (58) 

with the conventions of [Isj . 
The resulting action is 

S(P*,P) =£{E (T h [(V t - hb) + v] hIa (3 l2 
+ E JIM ( 

hhhli IM ) 

(59) 

where also for the 9j symbol we adopt the conventions 
of H|, and 

^iih = 7^ tr ( F h F i 2 e ) ~ 9h ^ 



WUisU = -2-90 + $>!« n tit(2^ + I)] 1 



/2 




(PhPiJiM = E (h,M 3 ,I 4 ,M 4 \I,M) p l3 M 3 (3 IiMi - 

M 3 ,Mi 

(60) 

The Hamiltonian is obtained 0] by omitting the 
time derivative and chemical potential terms, and re- 
placing the variables /?*,/? by corresponding creation- 
annihilation operators , b. These satisfy canonical com- 
mutation relations and should not be confused with the 



corresponding operators for the composites which are dis- 
tinguished by a hat 

Hb = E LO hM 1 i 2 M 2 b\ iMi b h M 2 + ^ ^ 



i 1 m 1 i 2 m 2 



I1I2I3I4 IM 



It is easy to check that, due to the symmetries of the 9j 
symbols, Hb is hermitian. 

From the interaction with external fields we get the 
fcrmion-boson mapping of other operators 



I 1 M 1 I 2 M 2 



M-2 



E (^) te( F hMiMFl M F l2M2 Fl u 



allz.Af 

xb l 1 M 1 b l 2 M 2 b hM 3 bl 4 M i - 



(62) 



Since the above Hamiltonian has been derived under the 
restriction n < < f2 in a single j-shell, we can assume 



-mim 2 — 6^mim2 3 



(63) 



so that the single boson energy matrix is diagonal 

wj lia = (2e-5/ i n)5/ lJa . (64) 

But the bosonic interactions couple all the bosons with 
angular momenta for which the 9j symbols do not van- 
ish, even if the corresponding potentials do vanish. It is 
however important to remember that the coefficients of 
the terms involving high angular momentum can change 
after the dependence of the index of nilpotency on the 
angular momentum is taken into account. 

B. Monopole pairing 

Let us consider the case of a pure monopole pairing 
interaction, namely gj = 0, 1 ^ 0. The exact spectrum is 

E niV — 2ne — go(ft + l)(n — v) + go(n — v) 2 + g Q {n — v)v. 

(65) 

We adopted definitions slightly different from the usual 
ones but more convenient, we think, in the present con- 
text. Here n is the total number of pairs, ie half the 
number of nucleons, and v is half the standard senior- 
ity, ie the number of pairs not coupled to zero angular 
momentum. It is then natural to relate the number of 
s-bosons n s to n and v according to 



n« = n — v. 



(66) 



In terms of n s , v, the eigenvalues can be written 
E ns +v,v = 2ne-g {fl+l)n s +g n 2 s +gon s {n-n s ). (67) 
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Setting bo = s and denoting by K angular momenta 
greater than zero we get from Eq. 1)61(1 the bosonized pair- 
ing Hamiltonian 



H 

+9o 



pairing 

K 1 ,K 2 ,K 



[2e — go(Cl + 1)] s + gos^ s s's 
(2Ki + 1)-?[{2K 2 + 1){2K 3 + 



x < K 2 0K 3 0\K 1 > 

IM 



(68) 



IM 



The expression of v\ is given in the next Section. The 
sector of zero seniority, n = n s , decouples and has the 
exact spectrum. The study of the seniority spectrum, 
in which the dependence of the index of nilpotency on 
the angular momentum might be important, is left for a 
future work. 



C. Monopole plus quadrupole pairing 

We consider now the case in which the nuclear Hamil- 
tonian contains a pairing plus quadrupole pairing inter- 
action. One might hope in a repetition of the pattern 
found in the bosonization of the pairing model, namely 
that the s — d bosons should decouple from the others, 
but this does not happen. We can write the Hamiltonian 
in the form 



H — H s -d + Hi 



Qi 



(69) 



(70) 



where Hq contains at least one boson with angular mo- 
mentum greater than 2 and i/ s -rf is the Hamiltonian 
truncated to the s — d subspace 

H s d — LdQ S + L)2 d) ■ d + VqS^ s s + w s d) ■ d + 

^vl{d) d))i ■ {dd) I +v 02 [(dU^)o,os s + h.c] 
i 

+60,2 [(rf t s t ) 2 ■ (dd) 2 + h.c] 

The parameters in the above equation are 

u>o = 2e — gofl 
uj 2 = 2e — g2^l 

J 3 

v = 2tt{ j j } g = g 


3 3 2 

50D { j j 2 } (2 9 2 - g ) 

2 2 1 

3 3 
w = 40ft <( j j 2 }g 2 

2 2 

3 J 2 

v 2 = 10ft { j j 2 } g 2 = V5g 2 




V02 



10V5ft 




go) 



y(3ff 2 -go)- 



(71) 



The last approximate equality holds for j 3> 1. Now 
in order to see if we have an at least approximate de- 
coupling of the s — d sector, we might proceed in the 
following way. After introducing the dependence of the 
index of nilpotency on the angular momentum, we should 
evaluate the energy of the s- and d-bosons first by assum- 
ing that the form factors of all the others vanish, then 
assuming nonvanishing only the form factor of the boson 
with angular momentum 4 and so on. If these energies 
get their minimum when all the other form factors van- 
ish, we have an exact decoupling. Otherwise we can have 
an approximate decoupling if these energies depend little 
on the inclusion of higher momentum bosons. We leave 
this investigation for a future work ^7j m which the de- 
termination of the form factors of the composites will be 
studied in a nucleon space including many shells. 



VI. THE PATH INTEGRAL FORMALISM: THE 
GOLDSTONE AND HIGGS FIELDS 

In this Section we outline the treatment of the effective 
action in the standard path integral formalism. This can 
be helpful in a numerical simulation. 

We must first determine the classical value of the s- 
boson field without breaking nucleon conservation. This 
can be done adopting the polar representation for this 
field 



Po = r exp I i — 
r 



r 



(72) 



(73) 



Boson number conservation is a consequence of the in- 
variance of the action under the transformation 

0o -> /3 exp(ia). (74) 

It leaves r, and therefore <r, invariant, while 

+ (75) 

We call ^3 an( i r the Goldstone and Higgs fields by 
analogy with the Godstone model, r is the classical field. 
The integration over a in the partition function extends 
from —r to oo, and one has to devise different approxi- 
mations depending on the value of r. All the other terms 
in the action must be expressed in terms of the Gold- 
stone and Higgs fields. For the time derivative terms, for 



9 



instance, we have 



-*rffVW0 + 5^-Vt/%f (76) 



where -FT refers to all the other bosons and VJ is the 
symmetric time derivative. The path integral must then 
be evaluated as a function of r and /xjv, and these pa- 
rameters are fixed by minimizing the free energy and by 
imposing condition J2J respectively. The evaluation of 
the path integral is performed by first determining the 
minimum of the action at constant fields, and then in- 
cluding the contribution of quantum fluctuations. But 
at the end we must check that 2n = n^. 



VII. SUMMARY AND OUTLOOK 

We have presented a new method of bosonization in 
which we restrict the partition function of the nucleus to 
nucleon composites. We get in this way the euclidean ac- 
tion of the effective bosons in closed form respecting all 
nucleon symmetries, in particular nucleon number con- 
servation. Indeed the presence of a large number of nu- 
cleons in the ground state is accounted for by a subtrac- 
tion which does not violate nucleon number conservation. 
The only approximation made concerns the replacement 
of the identity operator in the space of the composites by 
the operator V . 

The nuclear dynamics can be studied by the methods 
of path integrals, including numerical simulations which 
now are not affected by the sign problem, or in the more 
usual Hamiltonian formalism. Since the effective boson 
action does not have a canonical form, it has been nec- 
essary to devise an appropriate procedure to derive the 
boson Hamiltonian. This can be put in a form useful for 
practical applications, however, only by performing an 
expansion in the inverse of the index of nilpotency of the 
composites. 

The formalism is consistent only if all the composite 
bosons involved have a high index of nilpotency. Obvi- 
ously bosons with high angular momentum do not satisfy 
such condition, as shown by the example of a boson with 
maximum angular momentum 



&2j-l,2.;-l 



^2 (jwu'm 2 |2j -l,2j- l)c„ 



(77) 



which has index of nilpotency 1. It is then necessary to 
ascertain that bosons with low index of nilpotency will 
decouple. This problem sometimes is not explicit in some 
approaches, because the boson space is truncated directly 
or as a consequence of a truncation of the nucleon space. 
In this connection it is important to observe that the 



restriction to a definite boson space, like in the Interact- 
ing Boson Model, does not require an exact decoupling. 
An approximate decoupling is sufficient, because states 
weakly coupled can easily be integrated out in the path 
integral before deriving the Hamiltonian. This funda- 
mental feature is left for a future investigation. 

Another important point is the inclusion of particle- 
hole terms in the nucleon-nucleon potential. This re- 
quires the introduction of fonons, which will be done by 
an appropriate extension of the present technique. 

It is perhaps worth while to emphasize that for a 
derivation of the Interacting Boson Model it is not at 
all necessary that phonons be unimportant. It is suffi- 
cient that they can be integrated out like pions in the 
derivation of the nucleon-nucleon potential. 

In a future paper we will discuss in detail the determi- 
nation of the form factors of the composites 01 ■ This 
can be done by adopting the natural parametrization 

(BjM) mum2 = ^2pj a2 (jim 1 j 2 m 2 \JM). (78) 

The energies of the bosons are functions of the matrices 
p J which can be determined by a variational calculation. 
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APPENDIX A: BASIC FORMULAE FOR 
BEREZIN INTEGRALS 

In this Section we report for the convenience of the 
reader some basic formulae for Berezin integrals that we 
need. Their definition for a single Grassmann variable is 



^7(07 + b) = a, 



(AI) 



the generalization to many variables being obvious. For 
a change of variables 

7<=7i(Y) (A2) 
in a multiple integral we have 

/ nw/w = ( det w) 1 / nw/(7')- (as) 

Notice the appearance of the inverse of the jacobian, con- 
trary to the case of ordinary variables. 

Gaussian integrals can be evaluated exactly, like for 
ordinary variables. There are two types of such integrals 

/ \{{dl* h d lh ) exp ]T l*M ljlo = det M (A4) 

h ij 

J]J(d7h)expY,l^nj = PiA ( A5 ) 
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where Pf A is called the pfaffian of A. The following 
algebraic identity holds 



(PfA) 2 = det A. 



(A6) 



APPENDIX B: INNER PRODUCTS OF 
COMPOSITE STATES 

Let us first consider the case of only one composite. 
To evaluate the inner product of coherent states we use 
the identity operator in the fermion Fock space 



d7*d 7 ( 7 |7) 1 |-y) <Tl 



(Bl) 



where the 7*, 7 are Grassmann variables and I7) coherent 
nuclcon states 



I7) = exp(-7c t )). 



(B2) 



We then have 



(ftl/3) = = / d 1 *d 1 exp(-7*7)(/3 1 | 7 )( 7 |/3). 



(B3) 

The matrix element (/3i| 7 ) can be evaluated using the 
defining property of coherent states 



c|t) = 7I7) 



with the result 



(/3 1 | 7 ) = exp( ^=^7^7 



(B4) 



(B5) 



Therefore (/3i|/3) becomes 

(^1/3) = j d 1 *d 1 E{ 1 *, ll PlP), (B6) 
where the function E is 

1 



£(7*, 7, 0*,J9) = exp( -7*7+ /3* 7j B 7 

2v 



By the change of variables 

7=7- -ir( B ) 7 
p 

the integral is factorized according to 
(ft |/3) = | d 7 'exp('^=7 / /3 J Bt7' 
x dj exp 



(B7) 



(B8) 



~7 (Vn(/9i?t) *+ * ftflj ^39) 



The factors are of the form lA5l so that finally 



(/3i|/?) = 



det ( -^I3B^ 



det (v^(/?Bt)-i_ 



det 



1 + -=p 8{B^B 



(BIO) 



It is perhaps worth while noticing that in the limit of 
infinite f2, assuming the structure function to satisfy the 
condition l|38[) we recover exactly the expressions valid 
for elementary bosons, in particular 



1 + ^ 



exp(/3 t */3), ft^oc. (Bll) 



We will further need the inner product 

(/?!(&+)") = <?„(/?*)", 

where 



(B12) 



C 



1 



n — 1 



(B13) 

In the general case of many composites the above equa- 
tions become 



(ft|jS) = [det (t + $1*)]*, 



(B14) 



5"° <9" 



where 



exp |itrln[l + J=(a; • B*)*?]} | x=0 (B15) 
x-B^ = ^ Xj B\. (B16) 



APPENDIX C: THE OPERATOR V 

In order to implement the assumption of Boson Dom- 
inance we need the projection operator in the space of 
the composites describing the physical degrees of free- 
dom which dominate the partition function. Guided by 
the comparison with elementary bosons we approximate 
this operator by V . This comparison suggests that the 
measure in V should be (/3|/3) _1 . But, not surprisingly, 
we find that the measure must depend on the number of 
nucleons, namely of bosons 



d/*(/3*,/3) = J- 1 W) ~\ 



(CI) 



where the Jacobian is given in Eq. 1)28(1. 

For this reason it is necessary to introduce two types 
of composite bosons. We distinguish by a "bar" the new 
boson operators and their structure matrices. These new 
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bosons are to be regarded as the "physical" ones, while 
those appearing in V are auxiliary operators. 

Let us see the action of V on physical states. Let us 
first consider the case in which there is only one physical 
composite with structure function satisfying the equation 



B^B = 7 2 1. 



Using the results of Subsection B we find 

n 2 



(C2) 



v\(Vr) = j-^nf) 



(Q - n){Q - n - 1) 



(C3) 

Setting r 2 = n, V is the identity on states of n bosons 
provided 



1 /Q-n + 1 
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In the case of many composites, using the inequality 



(C4) 



(C5) 



which follows from the condition @ we find again that 
V approximates the identity with an error of order 1/57 

nb\ ) n °---(A) ni ) = I (Pl) n °-(bb ni +0(l/n)). (C6) 



APPENDIX D: DERIVATION OF THE 
EFFECTIVE ACTION 

For the following manipulations we need the Hamilto- 
nian in antinormal form 

H = i tt (j l+hQ y ch T c \_Y j g K ^cF K c^Fy (Dl) 

K 

where the upper script T means "transposed" and h was 
given in Eci. (|17|) . Now we must evaluate the matrix el- 
ement (/3 t | exp(—TiJ)|/3t_i). To this end we expand to 
first order in r (which does not give any error in the final 
t — ► limit) and insert the operator V between annihi- 
lation and creation operators 



(f3 t \ exp(-rfl)|A-i) = exp [~^tr(h + h )r j {(3 t \V 
~ch T TVc^9kT \ cF K cV \ C tF* C t|/?t_i). (D2) 

k 

Using the identity in the nucleon Fock space we find 
(/3 t |e X p(-rF)|/3 4 -i)= / d 1 *d 1 E( 1 *, 1 ,ff,[3 t _ 1 ) 



x exp [ —-tr(h + ho)r — 7*/ir 7 



x exp 



^5^7-7*47* J (D3) 

where the function E(j* ,7, j3* , (3) is defined in l|B7|l . By 
means of the Hubbard-Stratonovich transformation we 
can make the exponents quadratic in the Grassmann vari- 
ables and evaluate the Berezin integral 



((3 t \exp(-TH)\l3 t -i) = / Y[ da *K da K exp(-a* - a) 
J K 

x exp j itrln 1 + (^>* t + ^ ^g^a* Kl F K }j 



x det R exp ( — -tr(/i + ho)r ) , 



(D4) 



where 



R = 1 + h t. 



(D5) 



Performing the integral over the auxiliary fields a K * , a K 
we get 



(/3 t |exp(-rJI)|/3 t _i) = / da* K da K exp(-a* - a) 

^ K 

• j^trln 1 + (^>* t + V9^a* Kl F K }j 



x exp < -t 



x det R exp ( — tr(/i + ho)r ) , 



(D6) 



The functional form of the composites partition function 
is 



Zc 



dp* dp 



2ni 



exp (S eS (/?*,/?)) (D7) 



where S e Q is given in Eq. Qlfijl. 
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